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Motivated by experiments on atomically smooth layers of LaTiOa, a Mott insulator, sandwiched 
between layers of SrTiOa, a band insulator, a simple model for such heterostructures is studied 
using quasi one-dimensional lattices and the Lanczos method. Taking both the local and long-range 
Coulomb interactions into account, and computing the layer dependent local density of states, a 
metallic state was found at the interface whose extent strongly depends on the dielectric constant 
of the material. We also observed that the antiferromagnetic correlations in the bulk Mott phase 
persist into the metallic region. Our conclusions are in excellent agreement with recently reported 
results for this model in the opposite limit of infinite dimensions, thus providing an alternative tool 
to study electronic reconstruction effects in heterostructures. 
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I. INTRODUCTION 

The physics of nanostructures comprising atomically 
smooth layers of correlated electron materials with 
abrupt interfaces is a rapidly growing field spurred by 
experimental advances. Apart from the technological 
promise of new devices, these systems are of funda- 
mental interest for the novel phases that could occur 
at the interfaces of bulk systems. Correlated electron 
systems in the bulk are characterized by strong interac- 
tions and often exhibit a variety of competing states in 
close proximity within their phase diagrams. Fabricating 
dissimilar layers of these materials it has become possi- 
ble to achieve the enhancement of properties compared 
to the bulk, such as the polarization increase observed 
in superlattices made up of ferroelectric BaTiOs inter- 
spersed with paraelectric CaTiOs/SrTiOs, as compared 
to pure BaTiOa with identical volume^. In other experi- 
ments, coexistence of ferromagnetism and singlet super- 
conductivity has been observed in artificially fabricated 
superconductor (S) - ferromagnet (F) structuresSiS,, as 
well as the transfer of charge from a ferromagnet to a 
superconductor* . 

In this work, we shall be concerned with the reorgani- 
zation of charge and the interfacial state located at the 
boundary between a Mott insulator (MI) and a band in- 
sulator (BI), in the context of the experiments performed 
by Ohtomo et a&. These experiments involve the fabri- 
cation of a system with an arbitrary number of layers 
n of LaTiOa, a MI, sandwiched in between m layers of 
SrTiOa, a BI, on either side. The transport properties 
of this heterostructure (HS) revealed metallic behavior 
at the interface between the two types of insulators. Pi- 
oneering theoretical studies of this system by Okamoto 
and Millis&iZ used the Hartree approximation to treat the 
long range Coulomb interaction and a Dynamical Mean 
Field Theory (DMFT) to treat the local Coulomb re- 
pulsion. DMFT incorporates dynamical correlations ex- 
actly in the limit of infinite coordination and reduces the 
problem of each infinite layer to that of a single site cou- 



pled to a self-consistent bath. Consequently, the entire 
problem is mapped onto a one dimensional chain of inde- 
pendent impurity problems, coupled to each other via a 
self-consistency condition incorporating the hopping be- 
tween layers, and a Hartree mean-field for the long range 
Coulomb interaction. 

In this paper, we employ an alternative approach by 
considering a simplified one-band model for a finite het- 
erostructure using a quasi one-dimensional geometry. 
This limit of one dimension is the complete opposite to 
that of infinite dimensions that has been pursued in pre- 
vious literature for this model*^. Our main goal is to 
verify the conclusions of previous investigations, and pro- 
vide alternatives to DMFT for the study of HS systems 
involving correlated electrons. We treat the local and 
non-local Coulomb interactions exactly using the Lanc- 
zos method. Since the lattice parameters for SrTiOa and 
LaTiOa differ by only ~1.5% we neglect effects of strain 
at the interface and, as in earlier work^'^, we focus on the 
relevant electronic degrees of freedom within the Ti-band 
to define the Hamiltonian. We obtain the charge profile 
through the HS, along with the position dependent local 
density of states (DOS), and find a very good agreement 
with previous results^. In general, we observe that the 
strongest variation of the density moving from a MI to a 
metallic region, and then onward to a BI, occurs over a 
span of 3 to 4 unit cells. Furthermore, the antiferromag- 
netic (AF) correlations in the bulk Mott phase persist in 
the metallic interface. Similar results for the charge pro- 
file are found by extending our system to a 2-leg ladder. 



II. MODEL HAMILTONIAN 

Consider a HS with n layers of LaTiOa sandwiched in 
between m layers of SrTiOs. A one dimensional analog 
of this system is shown in Fig. [ija). La and Sr carry a 
bare valence of -1-3 and +2 within LaTiOa and SrTiOa, 
respectively. The condition of charge neutrality requires 
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FIG. 1: (a) Circles denote Ti-sites surrounded by La in the 
central Mott region and Sr in the band insulator (b) Charge 
density across the 16-site heterostructure for U = 20t and e 
varying from 3 to 15. 



that substitution of Sr by La in n-layers adds n electrons 
per unit area of each layer to the system. The essen- 
tial dynamics is then governed by the interplay of the 
electrons hopping on the Ti-sites, with their mutual long 
range repulsion, and the attractive potential produced by 
the surplus positive charge of the La ions. 

The Hamiltonian for this {n + 2TO)-sized HS in one 
dimension can then be written as 



where 



H = Ht + Hr. 
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Here Ht and Hu represent the nearest neighbor hopping 
and onsite Coulomb repulsion through the HS on the 
Ti-sites. Note that U is present on all the sites of the 
HS and not just in the MI region. The first term in 
Hcoul denotes the long range Coulomb repulsion between 
the electrons, with e being the dielectric constant of the 
host lattice. The second term represents the long range 
attractive potential between electrons on the Ti-sites and 
the La ions located at i?^ = Rj -f l/2a, with a being the 
lattice constant. The dielectric constant e is a strong 
function of temperature and frequency and also there is 
considerable uncertainty in the literature on the values 
of U for the system. Hence, we investigate the properties 
of the above Hamiltonian over a range of parameters, 
8t < i7 < 20i and 3 < e < 15 at r = 0, using the 



Lanczos technique^ with open boundary conditions. We 
set t — 0.3ey and a — 4A for simplicity. 




FIG. 2: (a) Coherent part of the total density, (b) Effective 
single particle potential and (c) Double occupancy across the 
heterostructure for U = 20t and e varying from 3 to 15. 



III. RESULTS 

Firstly, we present in Fig. [Ifb) results for the charge 
profile through a 16-layer (i.e. 16-site) HS with U—2Qt, 
as a function of the dielectric constant e. The structure 
is made up of 6 La-ions in the center and corresponds to 
a net free charge of 6 electrons. We shift the coordinate 
axis such that the center of the HS is z = 0. For e = 3, 
the charge density at the center is very close to 1 but 
it rapidly drops over four sites to a very small number, 
where the system approaches a BI. The intermediate re- 
gion where the density lies between 1 and is expected to 
be metallic. The strongest variation of the density mov- 
ing from a MI to a BI via a metallic region occurs over 3-4 
unit cells. Increasing e reduces the confinement of elec- 
trons in the central region within the MI and increases 
the width of the crossover metallic region. 

Since a metal is characterized by the presence of co- 
herent quasiparticles, we calculate the coherent fraction^ 
of the total density ricoh along the HS (Fig. [DJa)). This 
quantity is defined as an integral over the spectral weight 
from Lo=Q up to the first point in the occupied band where 
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the spectrum vanishes. We find that ricoh has a sharp 
peak in the crossover region reaching ^^90% of the to- 
tal density, decreasing to zero in the MI and BI regions. 
Both the results for the charge profile and ricoh are in 
excellent agreement with previous DMFT studies^. The 
density profile can be better understood by constructing 
an effective one-body potential V{z) which is obtained 
by substituting the result of the density across the HS 
for a given U and e into Hcoui and setting Hcoui equal 
to Y.^n^V{z). Fig.lDJb) shows V{z) for U=m,20t and 
e=3, 6. It can be observed that the effective potential 
that confines the electrons to the MI region in the cen- 
ter is practically independent of U and only depends on 
e. Increasing e reduces the depth of the potential well 
and reduces the slope of the potential in the crossover 
metallic region. 
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FIG. 3: Position dependent local density of states along the 
16-site heterostructure for U = 20t and e = 3. 

In Fig. mjc), we plot the double occupancy D along 
the HS for several dielectric constants, at a fixed value 
of U =20t. Note that D does not always follow the same 
behavior as the density through the HS. For e=3, where 
a clear MI phase prevails in the central region, D is sup- 
pressed at the center (z—O) relative to the maximum 
which occurs in the metallic state at the interface. Fi- 
nally D goes to zero as the BI is approached. Increas- 
ing the dielectric constant the double occupancy shows 
no dip at the center relative to the crossover region and 
essentially follows the behavior of the density, implying 



that there is a wider metallic phase that now evolves into 
a BI with no Mott phase in the center. 
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FIG. 4: Dependence of the local density of states on the di- 
electric constant for fixed density, n ~ 0.5, U — 20t. Lowest 
panel shows the same result for a uniform Hubbard model 
with periodic boundary conditions. 

In Fig [31 we show the position dependent local DOS 
along the HS for U=2Qt and e—3. At the center, z—0, we 
have a MI with density ri=0.999 reflected by a Mott gap 
separating the lower and upper Hubbard bands. Moving 
to the adjacent site (^=1) the density is only slightly 
lower at n=0.994 but a low weight in-gap state appears 
at the chemical potential. This metallic phase extends 
up to site 4, beyond which the DOS reflects a BI with 
all the spectral weight above a flnite gap approaching 
U. These results are in excellent agreement with those 
obtained earlier for this model using DMFT^. 

To clarify the nature of the metallic phase, in Fig.|4]we 
plot the local DOS for positions within the HS with com- 
parable densities (n«0.5) but different e, along with the 
spectral function for a quarter-filled uniform Hubbard 
chain of 16 sites under periodic boundary conditions. 
Both the partially occupied (lower Hubbard) and unoc- 
cupied (upper Hubbard) bands are significantly widened 
with the reduction of the dielectric constant and they de- 
velop long tails (although with low weight) as compared 
to the case with purely local interactions. For the low- 
est value e=3 most of the weight in the lower Hubbard 
band appears in a narrow feature arising from the onsite 
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Hamiltonian where the dielectric constant that governs 
the attractive interaction of the ions, Ca, differs from the 
dielectric constant for the electron-electron repulsion. 
In principle, there is no obvious reason why Ca — ^r- Our 
previous results simply followed Ref. 6 in these regards, 
but it is important to analyse the case 7^ e^- The mod- 
ified Coulomb term in the Hamiltonian then becomes, 



Hr 



ri+2m 9 



n-l-2m m+n 

E E 



The density profile through the HS for fixed 




(4) 



FIG. 5: (a) Spin-Spin correlation function along the het- 
erostructure. Inset shows the non-magnetic outer layers (b) 
Density profile in the 2-leg ladder also shows a metallic in- 
terface (c) Spin-Spin correlation function for the 2-leg ladder 
along the heterostructure for U = 20t, e = 3. 



Coulomb repulsion, which is well separated from the tails 
arising from the long-range part of that interaction. 

Next, we present in Fig. [U^a) the spin-spin correla- 
tions, (SozSjz) as a function of distance j from the cen- 
tral spin at 0. Clearly, AF order is favored in the central 
region with density close to 1, as is expected in a MI. As 
we move toward the outer layers of the HS, where the 
density is significantly reduced from 1 (beyond site 3), 
there is no alternation of sign in the spin-spin correlation 
function signifying a non-magnetic state as we approach 
the BI. However, in the intervening metallic interface we 
still have a clear signature of AF correlations as nearest- 
neighbor correlations (not shown) indicate. To check the 
robustness of our results we extend the system perpen- 
dicular to the interfacial direction by considering a 2-leg 
ladder (a 2x9 system with each leg containing 2 Sr ions 
on either side of 4 La ions) with U=20t and e=3. The 
density profile shown in Fig.[5fb) reveals a metallic region 
in between the MI and BI, as observed before for the ID 
case. In addition, from the spin-spin correlations along 
the ladder shown in Fig. [5Kc) we find AF correlations in 
the central MI. Nearest neighbor AF correlations persist 
in the metallic phase (not shown)^. 

Lastly, we consider a generalization of our model 



FIG. 6: Density profile through the HS as a function of ta for 
fixed — 6 and U = 20t. 

and varying is shown in Fig. [S] For < we see 
an increase in the density as we move away from the 
central Mott region, so that the density exceeds 1, and 
then drops across the metallic interface to approach zero 
in the BI region. When Ca > Cr we see a suppression of 
density below 1 in the central MI and a metallic region 
which increases in width at the interface. 



IV. SUMMARY AND CONCLUSIONS 

To summarize, motivated by experiments on 
LaTiOs/SrTiOa interfaces^, we have presented an 
exact diagonalization study of a simple one-band model 
for a correlated electron layered heterostructure in quasi 
one-dimensional lattices, where a MI is sandwiched 
in between a BI. Considering exactly the long range 
attractive Coulomb potential of the La ions and the re- 
pulsion between the electrons, we find that the interface 
supports a metallic state, in excellent agreement with 
previous studies of this model using DMFT and the 
Hartree techniques''^. 

The strongest variation of the density across the in- 
terface occurs over a region of 3-4 unit cells, whereas 
the actual width of the metallic region depends strongly 
on the dielectric constant of the material. The charge 
density depends only weakly on the local Coulomb in- 
teraction. Since the phenomena of interest occur over a 
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very short spatial range at the interface, we beheve that 
our results are valid even for HS much larger in the z- 
direction than the sizes we considered here. We further 
check the robustness of our conclusions by extending the 
system to a ladder configuration, finding similar results 
for the charge profile. The central region with density 
close to 1 shows AF correlations characteristic of a MI 
which persist into the metallic phase. 

The excellent agreement with DMFT doesn't imply 
that the self-energy for our system is relatively momen- 
tum independent, but that the observables such as charge 
density along the HS are unlikely to be strongly affected 
by momentum dependence of the self-energy. Cluster 
generalizations of DMFT such as cellular DMFT are 



likely to provide a better understanding of this issue. 
Important directions for future work using the method 
discussed here could involve realistic multiband models 
for other correlated electron phases including ferromag- 
netic and superconducting regimes. Also, the analysis 
of chain or ladder systems using density matrix renor- 
malization group or cluster generalizations of DMFT are 
possible, although the proper consideration of long-range 
Coulomb interactions are likely to pose a significant chal- 
lenge to these methods. 
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